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Abstract 
T^ ■ We recalculate the spectrum of radiation emitted from a collapsing star. We consider the simple model of a spherical 

f^ , star consisting of pressure-free dust and we derive the thermal spectrum via a systematic asymptotic expansion of 
^^ the complete Bogolubov amplitude. Inconsistencies in previous derivations are pointed out. 

(N 

More than a quarter of a century has passed since Hawking's remarkable suggestion [1] that a star collapsing to a 

00 black hole gives rise to radiation emission at a steady rate characterized by the black body spectrum. That particle 

^~~^ , creation takes place is not in itself surprising. Let us take the case of a spherical star with radius initially larger than 

its Schwarzschild radius that eventually collapses contracting to a point (according to classical gravitation). If we 

consider a quantized scalar photon field then the quantum spaces of the in states (before the initiation of collapse) and 

of the out states (after collapse is completed) are certainly different. Hence particle creation clearly takes place and 

r^^ is determined by the Bogolubov a{Lu, uj') and /3(w, w') amplitudes. The fact that the spectrum is that of a black body 

f^ is indeed noticeable, and ties up with the somewhat earlier results on black hole thermodynamics; see the textbook 

^O ' by Misner, Thorne and Wheeler [2] for a detailed treatment. Shortly after Hawking's work on black hole collapse two 

^P , classic papers by Fulling and Davies [3] , [4] appeared. The authors demonstrated an illuminating analogy, physical as 

well as mathematical, between gravitational collapse and the seemingly rather different problem of a perfect mirror 

starting from rest and accelerating for an infinite time; see the textbook by Birrell and Davies [5] for a comprehensive 

O I review of particle production by both mirrors and black holes. We shall see that the Fulling and Davies effect provides 

^^[ a very good guide as to what happens in gravitational collapse. At the root of both effects lies the fact that there is a 

?H . constant energy flux at late times as shown in both [1] and [4]. This does signal the presence of a thermal spectrum. 

PfJ' It also provides a physical interpretation of the logarithmic divergence that results upon integration of (28) over u;' 

^ in order to obtain the spectrum nluj) (the divergence will presumably be cured by a more complete theory). From 

the mathematical point of view Hawking's result is due to the special behaviour of the photon modes near the stellar 

surface just before the horizon is formed and his derivation makes heavy use of the asymptotic (i.e. near the horizon) 

j^ form of the modes. The singular behaviour of the modes in this regime has in turn given rise to various statements in 

" ■ the literature that are not strictly correct. For example there have been references to "parts of a{uj,uj') and (3{lo,lo') 

that relate to the steady-state regime at late times" ( [7] p. 327). However the Bogolubov amplitudes are global 

constructs and the distinction between early and late times does not make strict sense. The special role of late times 

does of course manifest itself in the calculation of local quantities like (7)xiy). The point of view that we are adopting 

is the one dictated by quantum-mechanical orthodoxy, namely that particle production is a global process ( [1] p. 

216) and that once initial and final states are specified (in the present case the initial state being the in vacuum) then 

the full standard expression for j3{uj,uj') is able to deliver the answer. Of course at some point in the calculation of 

(3{u,Lo') the special role of the horizon will show up but this must be taken care of by the mathematics without any 

ad hoc statements. In [6] we raised a criticism in connection to the calculation of the black body spectrum in [4] and 

the purpose of this note is to show that the points arising in [6] are also valid in the case of gravitational collapse. As 

in the case of the Fulling and Davies effect we shall perform a systematic asymptotic expansion of the amplitude for 

large w' and we shall show that as far as this particular calculation is concerned we can free ourselves from distinctions 

between early and late times, or transient versus steady state radiation. For an expanded version which examines in 

detail the points raised in this letter and which also to some extent duplicates results already found in the literature 

see [8]. Note: in what follows jh = G = c= 1. 

Our point is best illustrated in a concrete case and for the sake of simplicity we examine a collapsing spherical 
star consisting of pressure-free dust ( [9] chapter 11, sections 8 and 9). For the description of collapse one uses 
comoving coordinates (i, r, 9, ip) for the region inside the star and Schwarzschild coordinates (i, f , 9, ip) for the outside 
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region. Since the collapse is pressure-free each dust particle falls freely and thus the time t stands for the proper time 
registered by the particle in question. The two sets of coordinates are smoothly matched at the stellar surface. The 
time t registered by a surface particle and the radius R are given parametrically in terms of the cycloidal variable ( 
[9] equations 1.9.25, 26). The Schwarzschild time t for the surface particle is given in terms of the cycloidal variable 
by a more complicated expression (see [10] chapter 3, equation 107). It will prove crucial that wc specify the initial 
conditions of the collapse. Assume that the dust is held stable until time t — and is then allowed to contract under 
the action of gravity. Let M be the mass of the star, a its initial radius, and define a quantity k (as in [9]) by 

2M = ka^ (1) 

We shall consider the case of initially dilute matter a >> M, which implies via (1) 

ka^ « 1 (2) 

Hence at early times spacetime is essentially flat and the metric is simply expressed in terms of advanced and retarded 
Eddington-Finkelstcin coordinates v ~ i + f,u — i — r 

ds^ == dudv - f'^dn (3) 

(generally it is the Regge- Wheeler coordinate f* that enters in the definition of v and u but because of (2) there 
is little difference between radial and Regge- Wheeler coordinates at early times). At late times we need the metric 
expressed in terms of Kruskal coordinates; for the relevant definitions see pages 20-21 of [11] and also chapter 31 of 

[2]: 

r ' ' E ' \ 2M ) ^ ' 

The presence of E in (4) rcfiects an appropriate translation in Schwarzschild spacetime t and its precise value depends 
on the parameters of the collapse. 

We turn to the Klein-Gordon equation satisfied by a scalar field in curved spacetime. In this we follow section 
8.1 of [5] and p. 395 of [12]. A scalar mode corresponding to angular momentum I and satisfying the Klein-Gordon 
equation is written in the form 

i,i^{yA^ry' (t>i{r)Yr{0,^) (5) 

As far as conceptual purposes and technical details are concerned it suffices to restrict ourselves to s waves. Following 
the reasoning in [12] we neglect the centrifugal barrier (present even for s waves). Then the Klein-Gordon equation 
reduces to 

a-ii-aPi)*-» («) 

(in what follows the subscript in is suppressed). Since tl)o{r = 0) must be finite it follows from (5) that 

(j){r = 0) = (7) 

In terms of the Eddington-Finkelstein coordinates equation (6) takes the form of the massless scalar equation (8) of 
[6]. This, together with boundary condition (7) on the trajectory of the centre of the star, reduces the problem to that 
of a perfect mirror following a trajectory of the type examined in [6]. In particular the in and out modes ipLo{u,v), 
(fiuj{u,v) take the form of (15) and (16) of [6] respectively: 
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ipuj{u, v) = — {e^Tp{-iujv) - cxp i-iujp{u))) (8) 



1 

(fco{u, v) = — (exp {-iLofiv)) - exp {-iuju)) (9) 



Recall that apart from the modes ipi^{u,v) we should include modes quj{u,v) that contain no outgoing component; 
instead they are confined inside the black hole (the letter q complies with the notation in Hawking [1], [13]). The modes 



quj{u,v) are undetectable by an outside observer, but they are needed to make the set (punquj complete. Correlations 
between the two sets is part of the much discussed information problem [14], [15]. (Recall the analogy with the 
FuUing-Davies model: the mirror trajectory in [6] is asymptotic to the null line w = In 2 and thus does not cut all the 
characteristics of the wave equation; it leaves out the ones with u > ln2. The ensuing lack of completeness has been 
pointed out by Nikishov and Ritus [16].) 

To proceed with the calculation we need information on the function u = f{v) appearing in the expression for 
(pu){u,v). It will transpire that as far as large uj' asymptotics is concerned information pertaining to late and early 
times is sufficient. Consider an incoming ray l{i) (figure 1 of [8]) with v close to vh which upon reflection at the 
centre turns to an outgoing ray l'{i) with Kruskal coordinate U close to U ~ Q. For a given U we determine the 
corresponding u via the second of (4), and we then connect u of l'{i) to v of l{i) via /. Recall (see p. 29 of [11] for 
proof) that U is an affine parameter for a null geodesic U =constant. By a standard argument (p. 126 of [11], pages 
208-9 of [1]) we obtain 

. = /(.).-^ln(^),.^«, (10) 

We turn to the early stage of the collapse (figure 3 of [8]). Combining the expressions for t and R in terms of the 
cycloidal parameter we obtain for t small 

i? = a(^l-— j (11) 

Thus the star's surface initially contracts according to non-relativistic kinematics with acceleration M /a^ (in accor- 
dance with Newton's law of gravity). Given (2) the initial gravitational field is sufficiently weak so the light rays in 
Schwarzschild coordinates are straight lines. Let us consider a,t — R diagram and take an incoming light ray q which 
has w = 0. The reflected ray q' has w = and takes time a to travel from the centre of the star to the surface. During 
this time the gravitational held is still weak. This may be seen by setting t = a in (11); then (2) implies that by the 
time the ray emerges from the star's surface R still is almost equal to a. Thus at early times the function u = f{v) 
satisfles 

u = f (v) = V , initially (12) 

The Bogolubov amplitude f3{u!,uj') has the standard form (see e.g. [5], equations (2.9) and (3.36)) 

/3(w,w') = -«/ dzip^>{z,0)—ip^{z,0) + i / dz f — </?^'(z, 0) j ifiuiz^O) (13) 

where we use z for r* to ease up on notation and make contact with [6]. The integration in (13) can be over any 
spacelike hypersurface. Since collapse starts at i = the choice t = for the hypcrsurfacc is convenient. The in 
modes evaluated at i = are given by the simple expression (8) (i.e. p{u) = u) 

1 
(Puj{u, v) — {e'xjp{—iu!z) — cxp (iuz)) 6 (z) (14) 

where the presence of 9 (z) emphasizes the fact that z is a radial coordinate. The (p modes are given by (9) with / 
depending on the history of the collapse. Relation (13) is rearranged in the form 

/3(w,w') = ^ / " dz je'"^'^ - e'''^'^ 9 (z) {u^'e-'^^f - ujfe-"^f] + ^"^ ~^ / dz |e*'^'^ - e^^'^'^ (z) e*"^ 

ATT\/ujijj' Jo ^ J ^TTVLULj' Jo ^ ^ 

(15) 

The first and second terms in the above relation will be denoted by (3i{uj,uj') and /32(w, w') respectively. Observe that 
the argument of / runs up to vh and this is reflected in the limits of the first integral in Pi{lo,lo'). The rays with 
w < do not affect the amplitude due to the presence of the 9 function. Note that /32(u;,u'') is of kinematic origin 
and totally independent of the collapse. It was pointed out in [6] that (a) this term is unaccountably missing from [4], 
(b) that the term is in fact instrumental in obtaining the thermal spectrum in the case of an accelerating mirror. It 
will be seen that the same situation holds here (thus strengthening the analogy between Hawking and FuUing-Davics 
effects). One can see that the presence of the term /32(w,w') is indeed necessary by imagining that the star is held 
stable by some external means and that the radius stays constant (and large). Then trivially /(z) ~ z, the upper 



limit in the first integral in (15) is infinity, and the two terms identically cancel each other (as they must since nothing 
is produced). 

We return to the calculation of the amplitude and recall that we are chasing the leading asymptotic w' ^ oo 
behaviour of the amplitude which causes the logarithmic divergence mentioned in the Introduction. As in [6] we get 

/32(C.,C.')^— ^ (16) 

To calculate f3i{Lo,Lo') we perform an integration by parts to obtain 

Hiiio.io') = -^\^ r dze-'^f^'^-'^'^ (17) 

27r V uj J„ 

Observe that this is exactly the same amplitude that appears in (2.10b) of Davies and Fulling [4]. The leading order 
contribution will be examined using the method in [6]. The uj' ^ oo asymptotics of the integral in (17) belongs to the 
standard class of problems examined, for example, in Chapter 6 of Bender and Orszag [17]. To bring the singularity 
in the integral to zero we make the change of variable 

z — vh — X (18) 

and rewrite Pi{uj,uj') in the form 



(3,iu;,u;') = -— J- da;e— ^f-'e^ ^ (19) 



27r \ Lij jQ 

where the function g{x) = f (vh — z) is defined in the range < a; < vh and has the properties that follow from (10), 
(12): 

g(x)~-4Mln(|),a;^0 (20) 

9{vh) = /(O) (21) 

We momentarily drop prefactors and isolate the integral 

!=('' rfxe-^"ff(^)e^"'"^ (22) 

To obtain the asymptotic behaviour of (22) for w' large we use the same technique as in [6] (the two integrals being 
essentially identical). We deform the integration path to a contour in the complex plane; see Bender and Orszag op 
cit; chapter 6 of [18]; also [19], p. 610 where a very similar contour is used in the study of the asymptotic expansion 
of the confluent hypergeometric. The deformed contour runs from up the imaginary axis till iT (we eventually take 
T ^ oo), then parallel to the real axis from iT to iT + v^, and then down again parallel to the imaginary axis from 
iT + Vh to Vh- The contribution of the segment parallel to the real axis vanishes exponentially in the limit T ^ oo. 
We thus get 

POO POQ 

I = i dse-'^'^e"-'^^'^"^ -i ^g^^^'(v„+^s)^-^^g(vH+^s) (23) 



In both integrations the dominant contribution comes from the region where s ~ 0. Using (21) the second integral 
(including the minus sign in front) takes the form 

_if,i^v„^-iuf(n) dse-'^' = -i (24) 

Jo ^ 

In the first integral we use the asymptotic form (20) valid for small x and write 

e^pi~iujg{is))=c^pliAMuj\nl-jj = (-1 = cxp (^--4Mwj (^-^ 



where we took the branch cut of the function x^'^ to run from zero along the negative x axis, wrote x^' 
exp {iijj (In a; + i2N-K)) and chose the branch A^ = 0. Thus 



/•OO /'OO / TT 



(25) 



We substitute (24) and (25) in (23) and then in (19) to get 

e-"^'"" / TT \ ,4,,,,r(l + i4Mw) e-*'^^(°) , , 

27rvw^ ^ ^ ' [uo') i2it\/ijOlo' 

Cohecting (16) and (26) we get 



e 



-'""" „,,. ^ '^ . ,., A Z.-.4M. r (1 + z4Ma>) e— /(O) , 1 



/3 (cj, w') = -I— ^= exp - -4Mw ^-*4M^ ^ ,■■ ' ^= + ^= (27) 

The day is saved by (12) which causes an exact cancehation of the last two terms. The first term on its own 
immediately leads to the black body spectrum. Taking its modulus, squaring, and using the property 



\Y(\^iy)f = Try/ sinh (Try) 



we get 



,2 _ 4M 
2Tra;' e^ 



|/3(w,w')r = 5 — T7 (28) 



It is instructive to compare with earlier work on the subject. One often starts (see e.g. [5] p. 108, [4] equation 
(2.10b)) with expression (17) thus unjustifiably disregarding /32((-^,(-^') (in this respect cf the remarks following (15) 
above). One then uses (20) (or its equivalent (10)) throughout the range of integration and not just asymptotically 
near the horizon where it is valid. One thus gets dropping constant prefactors 

r^H , „ . iw4M , 
Pi{io,Lo')^J dxi^-) e-- (29) 

That (29) is an inadequate approximation to the original integral may easily be seen by the fact that the use of (20) 
has changed the behaviour of the integrand at a; = v^ (the non-singular end). From a physical point of view the use 
of (10) at z = induces a discontinuity in the derivative ,f'{v) at t = (corresponding in the FuUing-Davies model to 
an infinite acceleration of the mirror). One then proceeds to rewrite (29) by rescaling uj'x — > x 

Since one is chasing the ultraviolet divergence one simply sets uj'vh ~ oo, changes variable p ~ ia and rotates in the 
complex plane to get (30) in the form 

/>C30 

/3i(a;,a;') « ^--4M ^Jy^^^M-x ^-^^^m / ^^g-^^^- (31) 

Jo 

This reasoning is implicit in Hawking's statement ( [1], p. 209, lines preceding (2.19)) on the Fourier transform of 
the amplitude. Note that setting lj'vh = oo certainly does not amount to a systematic expansion in (w') . The a 
integration yields r(l + iuo) and one thus obtains the form for the j3i amplitude leading to the black body spectrum. 
On the other hand the step from (30) to (31) is again questionable. Integral (29) can be performed exactly in terms of 
the confluent hypergeometric function and the asymptotic estimate for large ui' may be examined afterwards. Indeed 
let us rescale the variable in (29) x -^ x/vh and rewrite 

l3l{0J,Uj') W ^;^4M+l^-,^4M f ^^^^u.'vHX^^u.4M ^ (33) 

^0 

= „«-'4M+i^-,.^4M^^^j (1 j^ -^^j^^ 2 + iivAM, iuo'vu) 
iuj + 1 



where M is the confluent hypergeometric. We can now examine the asymptotic hmit of (32) for large lo' . The 
asymptotic limit of the confluent M{a^ b, i \z\) for large values of \z\ is given by item 13.5.1 of [20] (z = iuj'vh)- In the 
case h = a + 1 some simplifications occur and we get 

M (1 + iw, 2 + ioj, i \z\) ~ - (1 + iLu) e'l^l -^ + ir (2 + iu) -^-j^ (33) 

\z\ \z\ 

(other terms are down by higher powers of 1/ \z\). The second term of the above relation combined with the pref- 
actors in (32) does feature the F (1 + ioj) e^^" factor characteristic of the black body spectrum. The reason for the 
discrepancy between (31) and (33) lies in the fact that one should first evaluate the integral in terms of the confluent 
and then take the lj' ^t oo limit rather than take the limit first. The rotation in the complex plane stumbles upon 
the Stokes phenomenon for the confluent (different limits for |2;| ^ oo depending on argz)). In short the black body 
spectrum (31) is obtained via steps identical to the ones in [4] and the criticism raised in [6] with respect to the latter 
apply to the case of gravitational collapse too. In this respect note the remark in [4] (the statement in parentheses 
following (2.10)) where doubts are implicit as to the validity of the aforementioned mathematical steps. One possible 
reason for the confusion may arise from a mathematical analogy with the Unruh effect (see e.g. chapter 5 of the book 
by Wald [21] for an excellent review). Equation (30) above is essentially identical to Wald's (5.1.11). Note however 
that Wald's (5.1.5) is exact for all V because of the kinematics of the problem and that the upper limit in (5.1.11) 
is strictly infinity because this is how the coordinates are defined. On the other hand the time t — Q (when collapse 
starts or when the FuUing-Davies mirror starts moving) and which is essentially responsible for the appearance of 
/32(w,tj') (16) simply does not appear in the Unruh problem (where the observer is forever accelerated). 
Acknowledgment 

The author is greatly indebted to Professor S A Fulling for his help. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7 



[9; 
[lo; 

[11 
[12; 

[is; 

[14 

[15; 

[16; 
[17; 

[is; 

[19 
[20 
[21 



S. W. Hawking 1975 Commun. Math. Phys. 43, 199 

C. W. Misner, K. S. Thorne, and J. A. Wheeler 1973 Gravitation (Freeman and Co., New York) 

S. A. Fulling and P. C. W. Davies 1976 Proc Roy Soc (Lon) A 348, 393 

P. C. W. Davies and S. A. Fulling 1977 Proc Roy Soc (Lon) A 356, 237 

N. D. Birrell and P. C. W. Davies 1982 Quantum Fields in Curved Space (Cambridge University Press, Cambridge) 

A. Calogeracos 2002 J. Phys. A: Math. Gen. 35, 3435 

B. S. DeWitt 1975 Physics Reports 19, 295 

A. Calogeracos 2003 "A detailed recalculation of the spectrum of radiation emitted during gravitational collapse of a 

spherically symmetric star", gr-qc/0406074 

S. Weinberg 1972 Gravitation and Gosmology (J. Wiley, New York) 

S. Chandrasekhar 1983 The Mathematical Theory of Black Holes (Dxford University Press, Oxford) 

P. K. Townsend 1997 ''Black Holes", gr-qc/9707012 

R. Brout, S. Massar, R. Parentani, and Ph. Spiendel 1995 Physics Reports 260, 329 

S. W. Hawking 1976 Phys. Rev. D14, 2460 

J. Preskill 1992 "Do Black Holes Destroy Information?", hep-th/9209058 

A. Heifer 2003 Reports on Progress in Physics 66, 943 

Nikishov A. L and Ritus V. I. 1995 Zh. Eksp. Teor. Fiz. 108 1121 

C. M. Bender and S. A. Orszag 1978 Advanced Mathematical Methods for Scientists and Engineers (IVlcGraw-Hill, New 
York) 

M. J. Ablowitz and A. S. Fokas 1997 Complex Variables: Introduction and Applications (Cambridge University Press, 

Cambridge) 

P. M. Morse and H. Feshbach 1953 Methods of Theoretical Physics (McGraw-Hill, New York) 

M. Abramowitz and I. A. Stegun 1972 Handbook of mathematical functions (Dover, New York) 

R. M. Wald 1994 Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics (University of Chicago 

Press, Chicago) 



